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The direct correspondence between Co band ferromagnetism and structural parameters is investi¬ 
gated in the pnictide oxides RCoPO for different rare-earth ions (R = La, Pr, Nd, Sm) by means of 
muon-spin spectroscopy and ab-initio calculations, complementing our results published previously 
[G. Prando et al, Phys. Rev. B 87, 064401 (2013)]. Both the transition temperature to the ferro¬ 
magnetic phase T c and the volume of the crystallographic unit cell V are found to be conveniently 
tuned by the R ionic radius and/or external pressure. A linear correlation between T c and V is 
reported and ab-initio calculations unambiguously demonstrate a full equivalence of chemical and 
external pressures. As such, R ions are shown to be influencing the ferromagnetic phase only via 
the induced structural shrinkage without involving any active role from the electronic / degrees of 
freedom, which are only giving a sizeable magnetic contribution at much lower temperatures. 

PACS numbers: 61.50.Ks, 71.15.Mb, 75.50.Cc, 76.75.-H 


I. INTRODUCTION 

In the last decade, quaternary RMX O oxides with 
tetragonal ZrCuSiAs structure (also generally referred to 
as 1111 compounds) have been extensively investigated 
on both experimental and theoretical levels. Here, one 
has an alternating stacked sequence of RO and MX lay¬ 
ers along the c axis and, within each layer, R 3+ (A 3- ) 
ions surround 0 2 ~ (M 2+ ) with a local tetrahedral crys¬ 
talline arrangement. Depending on the actual choice 
of the pnictide element X in combination with a tran¬ 
sition metal M, different electronic properties can be 
achieved^ ranging from high-T bulk superconductivity 
(with T values up to 55 K for M = Fe and X = 
to ferromagnetism with rare Hondo features (typically 
for M = Ru and X = P)L£ or with large - even colos¬ 
sal - magnetoresistance (M = Co, Mn and X = As, 
P),2r— However, the complex multi-orbital nature of 
the fermiology ; 16 ! 17 together with the open issue about 
the importance of electronic correlations ; 18 ’ 19 makes the 
overall understanding of these materials still not com¬ 
plete and currently highly debated. 

The role of the rare-earth element R is subtle as well. 
In the case of PFeAsO superconductors, it is well known 
that R strongly influences the optimal-doping T c value, 
which is more than doubled from ~ 25 K in the case 
of R = La to ~ 55 K for R = Sm. While it is ex¬ 
cluded that this enhancement is associated with 4/ or¬ 
bitals and with the localized magnetic moment at the 
R site^ no theory can still fully account for the ex¬ 
perimental observations. The structural origin (i. e., 
chemical pressure) seems to be a more adequate expla¬ 
nation for the observed effects^ however in this case 


it would be naively expected that external pressure in 
optimally-doped systems could enhance T as well, which 
is not the case experimentally! 22 ' 23 Chemical dilutions 
and quenched disorder are expected indeed to crucially 
drive the superconducting properties as well*2i invalidat¬ 
ing this way a full analogy between the effects of external 
and chemical pressures on T c 

In this respect, the impact of structural effects on 
the electronic ground state of the system seems to be 
much easier to control in the absence of intentionally- 
introduced chemical disorder. This was shown to be 
the case for undoped PFeAsO materials, where exter¬ 
nal pressure (P) weakens the long-range ordered (LRO) 
antiferromagnetism (AFM) similarly to what is achieved 
by means of low values of electron doping! 25 ' 26 At the 
same time, we have previously shown that the Co-based 
ferromagnetic (FM) compounds PCoPO are an ideal 
playground to investigate the interplay of chemical and 
external pressures^ 7 There, we reported by means of 
muon-spin spectroscopy (/i SR) and density-functional 
theory (DFT) calculations that ferromagnetism is size- 
ably enhanced with increasing P for both LaCoPO and 
PrCoPO. 2 ' In this paper, we extend our previous inves¬ 
tigation to the case of NdCoPO and SmCoPO. As it is 
well known in the literature, the magnetic moment aris¬ 
ing from electronic / shells and localized on the Nd 3+ 
and Sm 3+ ions plays a key role in the overall magnetic 
properties of the compounds^ - — similarly to the case 
of their PCoAsO analogues ! 12 ' 13,31 — In particular, a 
FM-to-AFM transition is induced within the Co sublat¬ 
tice with decreasing temperature, followed at much lower 
temperatures by a magnetic LRO phase within the R 
sublattice^ Here, we mainly focus on the FM phase at 
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higher temperatures by showing a strong correlation of 
the magnetic and structural properties. The critical tem¬ 
perature T c proper of the FM phase is shown indeed to 
be linearly dependent on the volume of the unit crystal¬ 
lographic cell and, in turn, on the ionic radius of R ions 
over a wide experimental range. At the same time, T c lin¬ 
early depends on the external pressure as well for all the 
investigated compounds in a quantitatively-comparable 
fashion. Accordingly, also supported by ab-initio calcu¬ 
lations, we provide a detailed and unambiguous evidence 
for a direct correspondence between external pressure 
and the lattice shrinkage induced by R ions. As such, 
these are demonstrated to affect the main ferromagnetic 
instability only via the induced structural shrinkage with¬ 
out involving any active role from the electronic / degrees 
of freedom. 


II. TECHNICAL DETAILS 

Loose powders of LaCoPO, PrCoPO, NdCoPO and 
SmCoPO were grown via standard solid-state reactions^ 
and their structural properties were investigated at am¬ 
bient temperature (T) and ambient P by means of x-ray 
diffractometry (XRD). The samples crystallized in the 
tetragonal phase, space group P4/nmm, diffraction pat¬ 
terns and Rietveld refinements being reported in Refs. [13 
andHil For the aim of clarity, the lattice parameters a 
and c for all the materials are summarized in Tab. [T] to¬ 
gether with the ionic radius value r l for the corresponding 
rare-earth ion R in a valence state 3+ and with a coor¬ 
dination number 8 (see Ref. [38:). A clear correlation of 
r, with the unit cell volume V is observed. It should 
be stressed that the chemical shrinkage of the cell driven 
by R ions is a well-known effect for the 1111 family of 
compounds j 31 i 39 i 40 

Measurements of dc magnetometry were carried out 
on a commercial Physical Property Measurement System 
(PPMS) and on a superconducting quantum interference 
device (SQUID) magnetometer Magnetic Property Mea¬ 
surement System (MPMS) by Quantum Design. 

fx + SR measurement s 41 ! 42 were performed at the Sp,S 
muon source (Paul Scherrer Institut, Switzerland), in 
conditions of zero external magnetic field (ZF) and for 1.6 


TABLE I: Results of Rietveld refinements of x-ray powder 
diffraction patterns reported in Refs. [2?] and 0. Lattice pa¬ 
rameters a and c for the investigated compounds are reported 
together with the corresponding volume V of the tetragonal 
unit cell and the ionic radius r T for the R ions. 


Compound 

a (A) c (A) 

V (A 3 ) 

r i (A.) 

LaCoPO 

3.968(7) 8.368(3) 

131.754 

1.18 

PrCoPO 

3.921(5) 8.212(4) 

126.253 

1.14 

NdCoPO 

3.904(2) 8.182(2) 

124.704 

1.12 

SmCoPO 

3.877(4) 8.072(9) 

121.331 

1.09 


K < T < 300 K. The GPS and Dolly spectrometers (7rM3 
and 7 tE 1 beamlines, respectively) were employed to per¬ 
form low-background measurements at ambient pressure. 
These were needed as independent references for the sub¬ 
sequent experiments performed on the GPD spectrome¬ 
ter (/xEl beamline). Here, external pressures P < 24 
kbar were applied at ambient T by means of a double¬ 
wall piston-cylinder pressure cell (PC) made of MP35N 
alloy, the transmitting medium Daphne oil 7373 assur¬ 
ing nearly-hydrostatic P conditions in the experimental 
ranged ac susceptometry was employed in order to de¬ 
tect the shift of the superconducting critical temperature 
of a small In wire inside the PC at T ~ 3 K and to ac¬ 
curately quantify P in turn, fx SR data anal ysis was 
performed similarly to what discussed in Ref. |27|, and 
the same notation will be employed in the current pa¬ 
per for the aim of clearness. In particular, the general 
expression 


A t ( t ) — A 0 



e -^-W + ( l -a PC )G s T (t) 


(1) 


was employed in order to fit the time (t) dependence of 
the measured spin (de)polarization P T = A T /A a for the 
implanted muons (/x + ) at all the investigated T values. 
Here, A T is the so-called asymmetry function while A 0 is 
an experimental spectrometer-dependent parameter ac¬ 
counting for the maximum value of A corresponding to 
full spin polarization (~ 100 %). When performing ex¬ 
periments on the low-background GPS and Dolly spec¬ 
trometers, a PC ~ 0 arises from the fx + implanted into 
the sample holder, into the cryostat walls etc. At the 
same time, for measurements on the GPD spectrometer, 
a PC ~ 0.5 0.65 mainly accounts for the fraction of in¬ 

coming /i + stopping in the PC. The spin polarization is 
here affected by the nuclear magnetism of the MP35N 
alloy resulting in a damping governed by the Gaussian 
and Lorentzian parameters cr PC and A PC . The T depen¬ 
dence of these quantities was independently estimated 
in a dedicated set of measurements on the empty PC. 
It should be remarked that, when ferromagnetic or su¬ 
perconducting samples are loaded into the PC and when 
performing measurements in the presence of an external 
magnetic field, a stray field is induced outside the sam¬ 
ple space. This leads to a complicated modulation of 
the local field within the PC which can be accounted for 
by introducing a proper profile-function and generalizing 
Eq. (JT]) accordingly! 23 ' 44 

In the case magnetic materials are being investigated, 
the function 


G s T (t ) 


[1 -VjT)]e-^ + 
a ± {T)F{t)D ± (t) + J(T)D^t) 


( 2 ) 


is employed in order to describe the spin (de)polarization 
arising from the remaining fraction (1 — a PC ) of fi , that 
is, the particles implanted into the sample. V m (T) is 
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FIG. 1: (Color online) Main experimental results for NdCoPO (left panel) and SmCoPO (right-hand panel) at ambient P: T 
dependence of the internal magnetic field at the fi site from the calibration low-background ZF-^t SR measurements (black 
filled symbols) and of the field-cooled dc magnetization M at H = 10 Oe (red filled circles). M for SmCoPO is reported after 
the subtraction of a small linear term accounting for the contribution of magnetic extrinsic phases. For both samples, the 
dashed gray line is a best fit according to Eq. m with £ = 0.34 as fixed parameter. Insets: ZF-fi SR depolarization curves 
at representative T values (ambient P, low-background spectrometers). Continuous lines are best-fitting curves according to 
Eqs. 0 and @. 


the fraction of n + probing static local magnetic fields 
(i. e., the magnetic volume fraction of the investigated 
material, due to the random implantation of /i + ). For 
V m (T) = 0, only nuclear magnetic moments can cause a 
damping of the signal. Typically this signal is of Gaussian 
nature with characteristic rate cr N ~ 0.1 /ts _1 . Below 
the critical transition temperature T c to the ferromag¬ 
netic phase, the superscript _L (||) refers to qi + probing 
a local static magnetic field in a perpendicular (parallel) 
direction with respect to the initial /r spin polarization 
and, accordingly, one has [a _L (T) +all(T)] = V m (T) for 
the so-called “transversal” (a- 1 ) and “longitudinal” (a^) 
fractions. A coherent precession of fi + around the lo¬ 
cal field B ' generated by a LRO magnetic order can be 
observed in the transversal amplitude. This is typically 
described by the oscillating function F(t), the standard 
choice F(t) = cos + being in good agreement 

with current experimental data (statistical y 2 ~ 1 — 1 . 2 ) 
at all T values. Here, 7 = 27tx 135.54 MHz/T is the gyro- 
magnetic ratio for fj + , while typical phase \<j>\ ~ 20° - 30° 
are measured (experimentally, negative values are found 
- see also the discussion in Ref. 1271) . The damping func¬ 
tion D 1 - (t) reflects a static distribution of local magnetic 
fields while the longitudinal component typically probes 
dynamical spin-lattice-like relaxation processes resulting 
in an exponential damping function D^(t). 

First principles calculations were performed using 
the VASP packag o 45 : 46 within the generalized gradi¬ 
ent approximation (GGA) to DFT— and projected 
augmented-wave (PAW)^ pseudopotentials for both the 
non-magnetic and ferromagnetic structures. Full relax¬ 
ation of the cell shape and of the internal positions was 
achieved using 650 eV as energy cutoff and a (12, 12, 6 ) 
fc-point mesh within the Monkhorst-Pack scheme!! and 


making the ab-initio forces to vanish (up to 0.001 eV/A) 
at each fixed cell volume. 


III. EXPERIMENTAL RESULTS 

A. Summary of main results from Ref. l27i on 
LaCoPO and PrCoPO 

In a previous workjll we reported on a fi + SR investiga¬ 
tion of magnetism in LaCoPO and PrCoPO. A LRO-FM 
phase with critical temperature T c = 33.2 ± 1.0 K and 
T c = 48.0 ± 1.0 K develops for LaCoPO and PrCoPO, 
respectively. The T dependence of the local magnetic 
field at the fi + site B M is well described for T < T c by a 
power-law functional 

B lt (T) = B lt ( 0)(l-^y (3) 

where £ = 0.34 ± 0.01 for both LaCoPO and PrCoPO. 
In spite of the qualitatively identical behaviour, La 3+ 
is a non-magnetic ion at variance with Pr 3+ . Accord¬ 
ingly, the only contribution to B (T) must arise from 
the Co sublattice alone in both the cases and, moreover, 
the functional form in Eq. © must be intrinsically asso¬ 
ciated to the FM phase from Co. The observed behaviour 
is qualitatively preserved independently on the P value, 
while a clear dependence on P was detected for both 
B II ( 0) and T c . A detailed discussion was presented in 
Ref. |27| concerning the former quantity. Here, we want 
to concentrate on T C (P), which was found to linearly 
increase with P in both LaCoPO an PrCoPO ^ 7 
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kbar (open symbols). Black filled symbols are relative to low-background measurements (reproduced from Fig. [I]). Insets 
display an enlargement of the T region just below T c . For both samples, the dashed gray lines are best fits according to Eq. © 
with ( = 0.34 as fixed parameter. 


B. NdCoPO and SmCoPO at ambient pressure 

The current n + SR measurements on NdCoPO and Sm¬ 
CoPO clearly evidence a different scenario if compared to 
what was reported for LaCoPO and PrCoPO. In fact, the 
magnetic moments of Nd 3+ and Sm 3+ ions crucially dom¬ 
inate the magnetic properties at low T and the interplay 
of d and / electrons leads to progressive re-ordering ef¬ 
fects on the Co sublattice . 28i3 ° This is confirmed by our 
measurements of ZF-/r SR and dc magnetometry both 
for NdCoPO and SmCoPO at ambient P. 

As shown in the insets of Fig. [TJ well-defined sponta¬ 
neous coherent oscillations appear in the p SR depolar¬ 
ization curves below ~ 55 K and ~ 65 K for NdCoPO and 
SmCoPO, respectively, indicating the development of a 
magnetic LRO phase in both the samples. These values 
also correspond to the sudden increase of the dc mag¬ 
netization M, as reported in the main panels of Fig. |T] 
Just below this transition, B^ (T) can be well reproduced 
by the power-law expression reported in Eq. m, again 
with £ ~ 0.34. By assuming that £ = 0.34 (fixed pa¬ 
rameter), a fitting of B (T) data leads to the estimates 
T c = 55.5±1.0 K and T c = 65.4T1.0 K for NdCoPO and 
SmCoPO, respectively (see the dashed lines in the main 
panels of Fig. [T]). However, such trend for B^{T) is no 
longer obeyed below T N ~ 16 K and T N ~ 45 K for Nd¬ 
CoPO and SmCoPO, respectively, where clear anomalies 
in the M(T ) curves are also detected. With further de¬ 
creasing T below T n , M vanishes while B (T) diverges to 
positive values in the case of NdCoPO, while a negative 
divergence is rather observed for SmCoPO [it should be 
recalled that the absolute value of B (T) is experimen¬ 
tally accessed indeed]. 

These results can be understood as follows. As B^ (T) 
is well reproduced by Eq. ©) with £ = 0.34 for T N < T < 
T c , it is reasonable to associate the magnetic transition 
at T c with a FM phase originating from Co similarly to 


what is achieved in LaCoPO and PrCoPO. In correspon¬ 
dence to T n another transition on the Co sublattice devel¬ 
ops with globally AFM correlations, as indicated by the 
gradual vanishing of M, possibly favoured by the mag¬ 
netic moments of the R ions. The positive and negative 
divergences of B^(T) in NdCoPO and SmCoPO, respec¬ 
tively, may be originating from several reasons among 
which, e. g., different microscopic magnetic configura¬ 
tions or statistical occupancies of the two /i crystallo¬ 
graphic sites in the two materials; 37 ’ 50 These issues are 
beyond the scopes of this paper and we will concentrate 
only on T c from now on. 

C. NdCoPO and SmCoPO for non-zero external 
pressure 

The impact of the crystallographic structure on mag¬ 
netism in NdCoPO and SmCoPO has been further inves¬ 
tigated through the application of an external pressure. 
As main result (see Fig. [2|, no qualitative changes are 
induced by P in the T dependence of B^ whose shape is 
preserved identical to what is reported in Fig. [T] How¬ 
ever, on a quantitative level, T c is enhanced by P in a 
linear fashion for both NdCoPO and SmCoPO, as en¬ 
lightened in the insets of Fig. [2] Overall, the situation is 
extremely reminiscent of what was observed in LaCoPO 
and PrCoPO, where P induces a linear increase of T 0 
as well. — Similarly to those compounds and as shown in 
the insets of Fig. [2j the shape of B (T) is not affected at 
all by P just below T c . Accordingly, Eq. 0 can still be 
used for T > T N as a fitting function with fixed parameter 
C = 0.34 in order to estimate T C (P) (see Fig.[5]later on). 
At the same time, T N is clearly suppressed upon increas¬ 
ing P for NdCoPO. In particular, T N = (16.25 ± 1.25) K 
at ambient P shifts to T N = (11.7 ± 0.7) K at P = 23.8 
kbar. The simultaneous T c enhancement and T N sup¬ 
pression is reminiscent of the effect of increasing mag- 
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netic field both in NdCoPO^ and in NdCoAsOJ^ On 
the other hand T N seems to be mostly unaffected by P 
for SmCoPO however, in this case the density of ex¬ 
perimental points does not allow to make any definite 
statement. We will no longer come back to the issue of 
the T n (P) dependence and only focus on T 0 (P) from now 
on. 


IV. DISCUSSION 

For itinerant ferromagnets, the magnetic instability is 
generally considered to be marked by the Stoner fac¬ 
tor S = I xc ■ D{E f ) becoming larger than 1. Here, I xc 
is the exchange integral, measuring the strength of the 
exchange interaction sensitive to the localization of the 
wave functions, while D(E F ) represents the density of 
states at the Fermi level, whose main contributions come 
from Co d-like components. An estimate of D(E F ) can 
be achieved on an ab-initio grounds at the experimental 
lattice constants of PCoPO at ambient pressure, as dis¬ 
played by the red squares in the main panel of Fig. [3] 
Accordingly, S can be computed as well by taking a 
value I xc = 0.99 for cobalt^i As a result, we find that 
the Stoner condition for FM instability is always safely 
satisfied for all the considered PCoPO compounds. In¬ 
deed, the FM state has a lower energy than the non¬ 
magnetic one, the total energy difference being of the or¬ 
der of 50 — 60 meV/cell for every compound. The green 
diamonds in the main panel of Fig. [3] are D(E F ) val¬ 
ues calculated for LaCoPO under external pressure and, 
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FIG. 3: (Color online) Main panel: density of states at the 
Fermi level D ( E F ) for the investigated PCoPO compounds as 
a function of the experimental value of the unit cell volume V 
(open squares, relative to lower x-axis). The D(E F ) values 
calculated for the LaCoPO compound are also reported at 
different volumes corresponding to pressure values indicated 
on the upper x-axis (open diamonds). Inset: dependence of 
the calculated unit cell volume in the FM state V calc on the 
experimental V for the four investigated PCoPO compounds. 
The dashed line represents the condition V calc = V. 
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FIG. 4: (Color online) Main panel: dependence of the critical 
temperature T c on the measured volume V of the tetragonal 
unit cell at ambient P for PCoPO and for the isostructural 
PCoAsO (these latter data are taken from Ref. 137") . The 
orange dashed curve is a linear best fitting to PCoPO data 
while the gray curve is a guide to the eye. Inset: dependence 
of the critical temperature T c on the ionic radius r 1 of P ions 
- see Tab. [I] The dotted curve is a linear best-fittings to data. 


R = Sm 



remarkably, these values very closely mimic the trend 
of the different compounds at their respective equilib¬ 
rium volumes, showing that the main variation is essen¬ 
tially linked to the internal pressure effect induced by the 
smaller R. This provides evidence that the change in the 
electronic properties of PCoPO is merely due to struc¬ 
ture and that the / states associated with the P ions do 
play a likely passive role. As shown in the inset of Fig. [3l 
the calculated equilibrium volumes V calc in the ferromag¬ 
netic state for all the considered compounds reproduce 
very nicely the experimental trend of V as a function of 
P (see the dashed line in the inset) both qualitatively and 
quantitatively within the usual DFT-GGA accuracy. The 
calculated magnetic moment of Co atoms decreases with 
reducing the volume of the tetragonal cell from 0.52 /i B in 
LaCoPO to 0.48 /r B for SmCoPO (not shown). While this 
trend is in qualitative agreement with the experimental 
results for LaCoPO and PrCoPO^i the strong magnetic 
contribution from P ions makes an experimental check by 
means of dc magnetometry uncertain for both NdCoPO 
and SmCoPO. 

We have shown above that, due to the effect of mag¬ 
netic P ions, the actual T dependence of B (T) for Nd¬ 
CoPO and SmCoPO strongly departs below T N from the 
behaviour reported for LaCoPO and PrCoPO. However, 
just below T c , strong qualitative analogies among all the 
samples suggest a common origin of the FM phase. To 
make these arguments more quantitative, we report in 
the main panel of Fig.[4]the T c values for the four PCoPO 
samples as a function of the relative measured V values 
for the tetragonal unit cell at ambient T and P and, in 
the inset, as a function of the r, of the corresponding P 
ion (see Tab. U). Not only it is evident that the lattice 
shrinkage triggered by smaller ionic radii is inducing a 
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FIG. 5: (Color online) Summarizing results of T c vs. P for the 
four investigated PCoPO samples. The slope of the dashed 
curve is obtained by Eq. 0 - see text. Data for PrCoPO, 
NdCoPO and SmCoPO have been shifted horizontally by the 
sample-specific amounts P biaB as discussed in the text. 
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dramatic increase of T c (being even doubled by a full 
replacement of La with Sm), but a linear correlation is 
clearly observed as well between T c and both V and r 1 . 

According to some theoretical models) 52 ’ 53 D(E F ) is 
expected to be directly related to T c and, as discussed 
above, we found that it slightly increases with the atomic 
radius of R and/or with external pressure. Although the 
increase of D(E F ) will surely result in a positive contri¬ 
bution to T c , we find from calculations that the small 
actual variation in its absolute value with chemical or 
external pressures can not explain the large differences 
reported in Fig. [7] for T c . Moreover, it should be re¬ 
marked that the peculiar linear relation between struc¬ 
tural parameters and magnetism is not at all a general 
feature of Co in 1111 oxypnictides. In particular, the 
main panel of Fig. [5] also reports T„ vs. V data for the 
isostructural compounds PCoAsO. 3 ' Here, the linear re¬ 
lation among these two quantities ceases as if the internal 
pressure effects are saturated yielding to a constant T c 
when approaching the smaller size of Sm. Accordingly, 
no change of T c has been observed as a function of ex¬ 
ternal pressure as well in SmCoAsO.— The differences 
between PCoAsO and PCoPO compounds imply a non¬ 
trivial difference in the microscopic origin of the link be¬ 
tween structure and magnetism. This may likely involve 
different degrees of hybridization between the pnictogen 
and the transition metal and, in turn, quantitatively dif¬ 
ferent exchange paths. 

Experimental results for the actual dependence of T c 
on P in PCoPO are summarized in Fig. 0 It is ap¬ 
parent that pressure leads to a linear enhancement of 
T c with a quantitatively similar slope for all the con¬ 
sidered PCoPO compounds. By considering the obser¬ 
vations above concerning the linear dependence of T c 
on V at ambient pressure, it is possible to discuss the 
close correspondence between chemical and external pres¬ 


sures more in detail. In particular, it is straightforward 
to speculate that external pressure linearly influences V 
which, in turn, linearly enhances T c . These arguments 
can be substantiated by computing the isothermal lattice 
compressibility /3 = — (1/P) x ( dV/dP ) for LaCoPO by 
means of first-principles calculations. The obtained val¬ 
ues are 7.52 x 10 -4 kbar -1 and 7.23 x 10 -4 kbar -1 in 
the paramagnetic and FM states, respectively. The as¬ 
sumption about the linear influence of both chemical and 
external pressures on V can now be exploited by writing 


1 dV 1 dV dT c 
VdP ~ ~VdT^~dP 


(4) 


and, accordingly, an estimate can be given for the partial 
derivative dT c /dP starting from the values computed 
above (the normalization factor 1/P is obtained from 
the ambient P value P = 131.754 A 3 ). This method 
is necessarily approximate, as P values at room T are 
involved together with T c , whose values are ~ 1 order 
of magnitude smaller. However, this uncertainty can be 
approximately solved in Eq. © by selecting an inter¬ 
mediate value for the computed compressibility, namely 
P ~ 7.4 x 10 -4 kbar -1 . This is admissible also in view 
of the tiny discrepancy (~ 4 %) among the two /3 values 
reported above. A line with slope given by the calculated 
dT c /dP is reported in Fig. [5] as well. One observes that 
experimental data are well-fitted by this line if the effect 
of chemical pressure introduced by R ions is mimicked 
by constant P bias offset values characteristic of each sam¬ 
ple. The experimental values P bias can be obtained from 
Fig. [5] by maximizing the agreement between experimen¬ 
tal data and the line. This linear correlation strongly 
supports the picture that R substitution acts as an ef¬ 
fective external isotropic pressure. A more stringent and 
physical argument to support this picture comes from 
the calculation of Pbias from first principles, with P bias 
defined as the calculated pressure needed to make La¬ 
CoPO match exactly the PCoPO equilibrium volumes. 
As shown in Tab. nn the excellent agreement between 
P 3 P and P 3 for all the considered samples is a clear 

bias bias x 

indication towards the full equivalence of chemical and 
external pressures in PCoPO. 


TABLE II: Experimental /P biaB J and computed /P biaB j val¬ 
ues of the bias external pressures mimicking the chemical pres¬ 
sure from R ions for all the considered compounds (see text). 


ex p 

Compound P biaB (kbar) 

theo 

P bias ( kbar ) 

LaCoPO 

- 

- 

PrCoPO 

55 

56.4 

NdCoPO 

78.5 

72.3 

SmCoPO 

107 

106.9 
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V. CONCLUSIONS 

In this paper, we have extended our previous p, + SR in¬ 
vestigation of itinerant ferromagnetism in LaCoPO and 
PrCoPO to the cases of NdCoPO and SmCoPO. The 
magnetic moment arising from electronic / shells and 
localized on the Nd 3+ and Sm 3+ ions plays a key role 
in the overall magnetic properties of the compounds at 
low temperatures. However, as long as only the ferro¬ 
magnetic phase at higher temperatures is considered, a 
strong correlation of the magnetic and structural proper¬ 
ties is reported and discussed. The critical temperature 
T c proper of the ferromagnetic phase is linearly depen¬ 
dent on the volume of the unit crystallographic cell and 
on the ionic radius of R ions over a wide experimental 
range. Its value can be even doubled by a full replace¬ 
ment of La with Sm. At the same time, T c linearly de¬ 
pends on the external pressure as well for all the investi¬ 
gated compounds in a quantitatively-comparable fashion. 
Accordingly, also supported by ab-initio calculations, we 
provide a detailed and unambiguous evidence for a direct 
correspondence between the structural shrinkage induced 
by external and chemical pressures. As such, the ferro¬ 
magnetic instability is driven in these compounds only 
by the crystallographic structure without involving any 
active role from the electronic / degrees of freedom. Cal¬ 


culations based on the variation in the density of states 
at the Fermi energy upon increasing pressures are not 
able to quantitatively reproduce the full enhancement 
T c , opening the way to further theoretical investigations 
of these materials. In this respect, we recall that 2D- 
3D crossover for spin-fluctuations have been put forward 
to explain previous experimental result o 30 ’ 31 and that a 
first-principle finite-temperature description of itinerant 
ferromagnets remains a big challenge in condensed mat¬ 
ter theory. 
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